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Abstract 

Purpose: To estimate tlie effects of lieterogeneity on tumour cell sensitivity to 
radiotherapy combined with radiosensitizing agents attributable to differences in 
expression levels of Epidermal Growth Factor Receptor (EGFr). 

Materials and methods: Differences in radiosensitivity are not limited to cells of 
different cancer histotypes but also occur within the same cancer, or appear during 
radiotherapy if radiosensitizing drugs are combined with ionizing radiation. A 
modified biologically effective dose (MBED), has been introduced to account for 
changes in radiosensitivity parameters (a and a/jS) rather than changes in dose/ 
fraction or total dose as normally done with standard biologically effective dose 
(BED). The MBED approach was applied to cases of EGFr over-expression and cases 
where EGFr inhibitors were combined with radiation. Representative examples in 
clinical practice were considered. 

Results: Assuming membrane EGFr over-expression corresponds to reduced 
radiosensitivity {aH = 0.15 Gy^ and aH/liH = ^-5 Gy) relative to normal radiosensitivity 
{a = 0.2 Gy^ and a//3= 10 Gy), an increased dose per fraction of 2.42 Gy was obtained 
through the application of MBED, which is equivalent to the effect of a reference 
schedule with 30 fractions of 2 Gy. An equivalent hypo-fractionated regime with a 
dose per fraction of 2.80 Gy is obtained if 25 fractions are set. Dose fractionations 
modulated according to drug pharmacokinetics are estimated for combined 
treatments with biological drugs. Soft and strong modulated equivalent 
hypo-fractionations result from subtraction of 5 or /(9 fractions, respectively. 

Conclusions: During this computational study, a new radiobiological tool has been 
introduced. The MBED allows the required dose per fraction to be estimated when 
tumour radiosensitivity is reduced because EGFr is over-expressed. If radiotherapy 
treatment is combined with EGFr inhibitors, MBED suggests new treatment strategies, 
with schedules modulated according to drug pharmacokinetics. 
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Background 

Recently, radiobiology has been transformed thanks to new knowledge concerning cel- 
lular activation processes in response to an external stimulus. This knowledge has led 
to the identification of promising new drug therapies called "targeted therapy" [1]. 

Epidermal growth factor receptor (EGFr) has emerged as a central molecular target 
for modulation during cancer therapy. The correlation between over-expression of 
EGFr and clinically aggressive malignant disease suggested that EGFr was a promising 
target for several epithelial tumours, which represent approximately two thirds of all 
human cancers. Furthermore, the favourable interaction profile for EGFr blocking 
agents combined with radiation has stimulated clinical trials in diverse anatomical sites 
including head and neck, colorectal region, pancreas and lung [2], where molecular in- 
hibition of EGFr signalling in combination with radiation represents a highly promising 
area [3,4]. 

Therefore, new radiobiology studies have focussed on identifying correlations be- 
tween radiosensitization and biological agents. However, these effects have not been 
fully integrated into current radiobiological models [5-8]. One such model commonly 
used in clinical practice, is the BED obtained from the LQ model [9], given by the fol- 
lowing equation (proliferation ignored): 



where a and ^ represent intrinsic and repair cell radiosensitivity, respectively, d repre- 
sents the dose per fraction and D is the total dose delivered during the radiation treat- 
ment. The BED is considered a "biological dose" delivered by a particular combination 
of dose per fraction and total dose to a given tissue, characterized by a given a/^ ratio, 
and is commonly used to equate or compare various fractionation schedules [10]. 

However, eq. (1) demonstrates that the same number of cells killed - the equivalent 
effect - could be obtained equating the BED not only for schedules with different num- 
bers of fractions and various doses per fraction, but also for schedules where the dose 
per fraction is increased if a reduction in radiosensitivity results (i.e. a or ^ is reduced). 

This could be applicable for subsets of cells that over-express EGFr, representing a 
source of heterogeneity closely connected with the repopulation rate and intrinsic 
radiosensitivity. However, the heterogeneous population of EGFr expression cannot be 
represented by a single equation of tumour control probability (TCP), as it is intrinsic- 
ally linked to a group of tumours with identical characteristics [11]. 

Furthermore, equations considering the radiation response that take into account dif- 
ferent compartments of sensitivity within tumours [12] or a Gaussian distribution of in- 
dividual radio sensitivities [13,14] cannot be used because various levels of 
radiosensitivity coexist in the tumors or in the statistical sample. 

Therefore, during this computational study, a new mathematical interpretation of 
radiosensitivity parameters that are normally used in standard radiobiological models 
(i.e. as functions of EGFr expression) is proposed using simple examples. 

The final aim of the current study is to provide an additional mathematical tool that 
can be used to carry out radiobiological analysis, taking into account the radioresis- 
tance effects due to EGFr over-expression and/or radiosensitization effects due to EGFr 
inhibitors when they are combined with radiation. 




(1) 



Pedicini et al. Theoretical Biology and Medical Modelling 201 2, 9:23 
http://www.tbionned.conn/content/9/1/23 



Page 3 of 1 2 



These examples are not intended to simulate a particular type of radiotherapy treat- 
ment, but are designed to demonstrate a general effect. 

Materials and methods 

During the current analysis two separate groups of patients with various levels of EGFr 
expression were considered. For each of the EGFr groups, various values for the para- 
meters a, ^, were considered. This approach allowed radiobiological analysis to be con- 
ducted in cases where differences in radiosensitivity occurred within the same tumour 
after combined treatments comprising radiation and radiosensitizing EGFr inhibitors 
[4,15,16]. In the latter case, various levels of radiosensitivity did not coexist, but they 
followed one another according to the concentration of radiosensitizing drug present 
during the radiotherapy session (Figure 1). 

Modified BED: Effects due to a change in EGFr expression levels EGFr expression has 
been assessed through intensity of staining (i.e., absent, minimal, moderate, or intense 
staining) in clinical practice [17]. During the present analysis, normal and high expres- 
sion levels of EGFr (i.e. below and above 50% staining) were distinguished. The sub- 
script H was added to indicate high EGFr expression. 

The BED for the EGFr group with high expression may be indicated as: 

BEDh = n'd(l + 



Here, because an and are lower than a and ^ (reduced radiosensitivity), the num- 
ber of cells killed with the same dose per fraction {d) and the number of fractions {n) 
were reduced with respect to standard radiosensitivity conditions. Therefore, the fol- 
lowing inequality arose: 

a BED > an 'BEDh 

To obtain the same effect with an equal number of fractions, a change of dose/frac- 
tion is necessary. We introduce the MBED: 

MBED = n-8- (1 + —^r^ (2) 

where the dose 8^ which refers to an and y^//, has the effect equivalent to d, which refers 
to a and ^, so that: 

a BED = an MBED (3) 

In eq. (3) the LHS provides a measure of treatment effect under standard conditions 
of radiosensitivity, while the RHS represents the same effect achieved under non stand- 
ard conditions of radiosensitivity. 




Figure 1 Schematic representation of radiosensitivity variability within a single tumour due to the 
presence of varying concentrations of radiosensitizer drugs (Light gray = high radiosensitivity, dark 
gray = low radiosensitivity). 
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From eq. (1), eq. (2) and eq. (3): 
a-d' ( 1 + ) = an-S' ( 1 + ^ 



and solving for 6 



Therefore, the MBED distinguishes between changes in biological effect due to irrep- 
arable and/or reparable damage variations, rather than changes due to dose/fraction or 
total dose variations. A reduction in radiosensitivity due to increased membrane EGFr 
expression [11,18] implies equivalence between treatments by increasing the dose per 
fraction with an equal number of fractions. 

Furthermore, to obtain isoeffective treatments with a different number of fractions m 
{m<n hypo-fractionation, m>n hyper-fractionation) from eq. (3), the following results: 



(5) 




Modified BED: Effects due to biological drugs 

Combined treatment comprising radiation and radiosensitizing EGFr inhibitor drugs 
requires the daily dose that achieves the same effect without drugs to be calculated. 
This will result in a calculation of the daily radiosensitivity conditions induced by the 
drug compared with standard radiosensitivity. 

On the basis of a preclinical assessment, we propose a method to estimate the daily 
radiosensitivity when radiotherapy treatment is combined with biological drugs. Subse- 
quently, the MBED method is applied to assess the changes required in terms of dose 
fractionation when such daily radiosensitivity is considered. 

During the first phase, survival curves obtained with various concentrations of a 
monoclonal antibody (mAb) EGFr inhibitor were selected from the literature [16,18]. 
From these curves, using a polynomial regression, the corresponding values of a and ^ 
were calculated (Figure 2(a)). However, the drug concentrations reported in these 




IGy 2Cv 3Gv 4Gy 5Gy MOll TUC WCCl THu Fli 



Figure 2 First phase to investigate the effects EGFr over-expression on radiosensitivity of Head and 
Neck cell lines. Data from literature [16,18,19] demonstrate the correlation between EGFr over-expression 
and reduced cellular radiosensitivity. This situation is indicated by an upward shift of the cell survival curve 
in the line over-expressing EGFr compared with normal EGFr expression. A polynomial regression allows 
radiosensitivity parameters corresponding to various surviving curves to be calculated. 

V ) 
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studies do not correspond to the effective drug concentrations used during the com- 
bined treatment with radiation every day of treatment (Figure 2(b)). 

Therefore, during the second phase, the daily in vivo concentration of EGFr inhibitor 
drug was calculated from its pharmacokinetic curve and drug dosage [20]. Referring to 
these daily concentrations, it is possible to interpolate plausible corresponding curves 
of survival fractions, obtaining the researched values of a and ^ using a new polynomial 
regression (Figure 3). 

Subsequently, assuming a daily in vivo radiosensitivity, eq. (3) with a variable con- 
centration of a radiosensitizing drug according to the weekly dosage can be written 
as follows: 

n^dia ^ ^-d) = n^'[Siiai ^ P^'Si) ^ + 65- {as +^5-^5)] 

where {n^ = m/5) represents the number of weeks of overall treatment and the 
numbers 1,2,. . .,S indicate the day of the week. In compact form, we can write: 



i=l 



Therefore, an equivalent fractionation with the same number of fractions is obtained 
using the following: 



1 ^ 



(6) 



From eq. (6), a solution with equal dose for each day is: 



(7) 



In addition, eq. (6) highlights the possibility of solutions with a dose adapted to the 
daily radiosensitivity. By equating the effect day to day during the week, for the ith day 
we obtain: 



d-{a + ^-d) = 8i 



-Mi, 




a ? = I6G) flOnXf) 
a $=IS G} (13 nM) 
$ = HGy(:2nM) 
l!Gy(}7nM) 
^=9G} (61 nM) 



Mon Tue Wed Thu Fri 



Figure 3 Second phase to estimate the effects on radiosensitivity of variable concentrations of 
mAb EGFr inhibitor in Head and Necic cell lines. Surviving fraction curves corresponding to tlie daily 
concentrations of mAb from pliarmacol<inetics curves [20] are obtained by interpolation. The following 
concentrations of EGFr mAb inhibitor are obtained: 100, 61, 37, 22 and 13 nM. The corresponding 
polynomial regression curves provide a/p> = 5, 9, 12, 14 and 15 Gy (with a = 0.2 Gy with respect to 
untreated cells with a/p= 16 Gy {a = 0.2 Gy). 
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therefore: 

^'-w^MMff) 

Eq. (8) leads to a modified fractionation modulated according to the pharmacokinetics of the 
drug combined with radiation. For a schedule with different numbers of fractions : 

This solution leads to a modulated hypo-fractionation if the number of weeks is less 
than the standard fractionation (vice versa for the hyper- fractionation), 

Eq. (7) and eq. (9) represent dose values that have the same effect. However, as 
in the drug is also absorbed by normal tissue cells, these cells will show increased 
radiosensitivity. Therefore, modulated dose fractionation with a reduced dose of ra- 
diation corresponding to higher radiosensitivity could lead to a reduction in harm- 
ful effects. 

This proposal could be verified through clinical trials. 
Results 

This section discusses results from representative examples occurring in clinical prac- 
tice. Schedules with the equivalent effect of 30 fractions of 2 Gy/fraction (assumed as a 
reference standard regime) were calculated. To analyze an increase in radiosensitivity, a 
change in a or yS, and consequently a change in a/^, has been assumed to simplif)^ the 
calculations without losing generality. 

For examples 3, 4 and 5, substantial changes in ^ alone has been adopted, assuming 
that data were obtained from the polynomial regressions of curves depicted in Figure 3. 

Of note, the unchanged a, /3 (without polynomial regression) and the fractionation 
schemes assumed in these examples are plausible but should not be considered as 
recommendations for real clinical situations. 



NozimIECFR (a/J}=10qy) 



2Gy 



(a) 



Radio*ensiiivity 



Hicli. 



ECFR Ov«i«xpi«<««d (a/^=7.50y) 



Dose 

3Cr- 
2Cy- 
ICy. 



ECFR Ovei«3ivz«ft«a (<f/fi=7.50y) 



Weekljr dose firactioitation 



Weekly dose fractionation 



Weekly dose fractionation 



Weekly dose fractionation 



Figure 4 Weekly dose/fraction as a function of radiosensitivity for modified fractionations with (a) 
same number of fractions as the reference fractionation {Example 7) and (b) hypo-fractionation with 
one weeic less than reference fractionation {Example 2). 
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Radiosensitvity 

High 

Medium 

Low 
Dose 



6 weeks 



(a) 



5 weeks 



o 3 



2^Gy 



(b) 



4 weeks 



o2 



.3 



(c) 



Mon Tuc Wen Thu Fri ' Mon Tuc Wen TTiu Fri ' Moit Tue Wen TTiu Yri ' 

Figure 5 Dose/fraction as a function of daily radiosensitivity for modulated fractionations with (a) 
same number of fractions as the reference schedule {Example 3: 6 weeks), (b) hypo-fractionation 
with one week less than the reference schedule {Example 4: 5 weeks) and hypo-fractionation with 
two weeks less than the reference schedule {Example 5: 4 weeks). The grey lines represent 
radiosensitivity corresponding to tlie pliarmacokinetics curves of absorption for tine EGFr mAb inliibitor. 
Abbreviations: (1 ) a//3 = 5 Gy; (2) a/p = 9 Gy; (3) a/^ = 1 2 Gy; (4) a//3 = 1 4 Gy; (5) a//3 = 1 5 Gy. 



Dose fractionations are presented for examples demonstrated in Figure 4 and Figure 5. 
These figures present the extent of dose for fraction as a function of weekly or daily 
radiosensitivity; tables 1 and 2 present numerical results. 



Example 1 

In this example a selection of patients that should be treated with the reference sched- 
ule (consisting of 30 fractions of 2 Gy/fraction on PTV) was assumed. Patients in the 
first subset (Gl) were considered to have normal EGFr expression on clonogenic 
tumour cells, with radiosensitivity corresponding to a = 0.2 Gy^, ^ = 0.02 Gy^ {a/^ = 10 
Gy), In addition, we considered a second subset of patients (G2) as presenting with 
EGFr cell membrane over-expression, resulting in a reduction of radiosensitivity with 
an =0.15 Gy\ Ph=0.02 Gy'^ {a/^n^ 7.5 Gy). 
Therefore, with respect to the reference schedule, the effect for the subset Gl would be: 

aSED = 0.2.30 2. (1 + ^ | = 14.4 

Owing to the reduction of the a component of irreparable damage, the same schedule 
used for group G2 will produce the following effect: 

a^.^^DH = 0.15.30.2. M+;^j =11.4 



with a noticeable reduction in the effect of overall treatment. 
Table 1 Numerical results for Examples 1 and 2 



EGFr expression 


a{Gy') 


i8(Gy-') 


amOy) 


dexAGy) 


dex2{Gy) 


Normal 


0.2 


0.02 


10 


2.00 


2.33 


Over-expressed 


0.15 


0.02 


7.5 


2.42 


2.80 



Abbreviation: d^xi and dex2 = doses from MBED for Example 1 and 2, respectively. 
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Table 2 Numerical results for Examples 3, 4 and 5 



Day 


a{Gy') 




a/PiGy) 


dexsiGy) 


dexAGy) 


dexsiGy) 


Monday 


0.2 


0.040 


5 


1.68 


1.94 


2.31 


Tuesday 


0.2 


0.022 


9 


1.86 


2.18 


2.62 


Wednesday 


0.2 


0.017 


12 


1.93 


2.27 


2.74 


Thursday 


0.2 


0.014 


14 


1.98 


2.32 


2.81 


Friday 


0.2 


0.013 


15 


1.99 


2.34 


2.84 



Abbreviation: dexs, dex4 and c/ex5 = doses from MBED for Example 3, 4 and 5, respectively. 



To produce the same therapeutic effect for patients in G2 as received by patients in 
group Gl (with the same number of fractions taken in the reference treatment), the dose 
per fraction should be increased by imposing condition (3). Then, from eq. (4), we obtain: 

6 = 2A2Gy 

To achieve the same effect on the PTV, 30 fractions of 2.42 Gy/fraction should be 
given to compensate for reduced radiosensitivity due to over-expression of membrane 
EGFr (Table 1 and Figure 4). 

The new schedule will be not equivalent in terms of toxicity to organs at risk (OAR). 
Therefore, the plan will require re evaluation of the harmful effects for OARs. In the 
opposing situation, that is for an increase of radiosensitivity in the clonogens of G2 
compared with Gl (owing to a radiosensitizing drug), one can adopt the same proced- 
ure. In such cases, the equivalent effect on the PTV, with the same number of fractions, 
will be reached by reducing the fraction dose. 

Example 2 

For the same subsets of patients used in Example 7, we analyzed a hypo-fractionated 
schedule that lasted for one week less for patients in G2, with the same effect as the 
standard schedule for patients in Gl. In the hypo-fractionation schedule, the number 
of fractions was m = S-(n^-l) = 5-5 = 25 fractions. 
Applying eq. (5) we obtain: 

8 = 2.80Gy 

therefore, the hypo-fractionated schedule for patients in G2 will be equivalent to the 
standard schedule for patients in Gl if 25 fractions of 2.80 Gy/fraction are given. If 
a/^-10 Gy and a normal radiosensitivity is assumed, we would obtain: 

0.2.30 2. 1 1 + ^ I = 0.2.25.^. ( 1 + -^V 

from which: 

d = 233Gy 

which would underestimate the dose required to achieve the same effect on the 
PTV (Table 1 and Figure 4). 
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Example 3 

In this example we refer to group G2 having substantial membrane EGFr over-expres- 
sion, with an =0.2 Gy'^ and Uh/^h^ 1^ Gy (similar estimated a/^ values are reported in 
the literature [21]), We compare the reference treatment with a combined treatment 
comprising radiation and biological drugs that produce an increase in radiosensitivity. 

In addition, we assume a weekly drug dosage with a pharmacokinetics curve showing 
maximum absorption during the first day of treatment [20]. The weekly radiosensitivity 
is assumed to be that described by the set of values reported in Table 2. 

The equivalent treatment with the same number of fractions is obtained using eq. (7). 
In this case, a constant dose for each day is obtained, equaling the global effect. 

6 = l.88Gy 

Subsequently, using eq. (8), a dose modulated according to the drug pharmacokinet- 
ics is obtained, equaling the effect for each day of treatment. Results are presented in 
Table 2 and Figure 5. 

Examples 4 and 5 

The equivalent global effect of the reference schedules could be obtained by subtracting 
one or two weeks of treatment from eq. (9), with a modulated soft hypo-fractionation 
(5 weeks) and with a modulated strong hypo-fractionation (4 weeks), respectively. The 
results are presented in Table 2 and Figure 5. 

Discussion 

During practical applications of radiobiological models, the main difficulty is to decide 
which parameter values should be included in individual calculations. It is important to 
clarify that population based estimates of the a/^ value represent averages, and that 
values are likely to vary between and within tumour types. It is clear that the assump- 
tion of a single value for a or a/^ is a simplification and this could have a considerable 
impact on the predictive use of BED when deciding on dose fractionation [22] . 

However, recent knowledge concerning molecular mechanisms allows new develop- 
ments to be explored and provides important information relating to the intrinsic 
radiosensitivity and fractionation sensitivity. Cell studies in vitro demonstrate that dif- 
ferences in radiosensitivity occur among cell lines derived from different types of 
tumours or from the same type of tumour, and during irradiation when combined 
treatments using radiation and radiosensitizing drugs are utilised [16,23-25]. 

These considerations may lead the way for new studies concerning evaluation of a 
and yS, in which cellular radiosensitivity is modified using known concentrations of 
radiosensitizing drugs, as described in Figure 4 and Figure 5. 

Therefore, the historical inability to distinguish among effects resulting in differences 
in radiosensitivity could be overcome through new knowledge concerning heterogen- 
eity [26,27]. These effects are well known from preclinical studies, and could be used to 
reduce uncertainties and investigated through clinical trials [28]. The ideal situation 
could be to use assay methods to allocate patients to various treatment schedules on 
the basis of individual measurements of tumour cell radiosensitivity (for example, due 
to varied expression of EGFr) or absorption of drugs. This approach is expected to be 
applied in the foreseeable future. 
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On the basis of these considerations, a new method to interpret BED expression, 
named MBED, was introduced during this computational study to take account of in- 
trinsic differences in radiosensitivity. 

The requirement to introduce MEED arises because radiosensitivity is usually consid- 
ered to be fixed for a cell type and constant during any radiation treatment. For this 
reason, a and ^ are considered fixed values with considerable uncertainty. Therefore, in 
the standard use of the BED, the hypothesis that one fractionation is equivalent to an- 
other underlies the assumption that the values of a and /3 are the same: to have the 
same effect - resulting in the same number of cells being killed - changing the dose 
per fraction, one must alter the number of fractions. 

Herein, it is argued that for various values of radiosensitivity, the same number of 
cells can be killed with the same number of fractions by varying the dose per fraction. 
This requires identification of prognostic parameters such as the over-expression of 
EGFr, which allows the radiosensitivity of the individual patient to be classified and the 
most appropriate radiation dose fractionation to be identified. 

The results of this study demonstrate that for a subset of patients presenting with 
EGFr cell membrane over-expression, resulting in reduced radiosensitivity with respect 
to a subset of patients with normal EGFr expression of clonogenic tumour cells, the 
dose per fraction should be increased to produce the same therapeutic effect with the 
same number of fractions taken in the reference treatment. 

When radiation is combined with a biological drug that produces an increase in radio- 
sensitivity, depending on the drug dosage, the equivalent treatment with the same num- 
ber of fractions is obtained by a dose of radiation modulated according to drug 
pharmacokinetics. 

The dose needs to be increased if the number of fractions is reduced. 

In the examples reported herein, the absorption of EGFr inhibitors was considered for 
cancer cells alone. In general, cells of normal tissues also absorb the drug. In particular, 
EGFr is over-expressed in skin cells. Therefore, the effect of increased radiosensitivity 
will affect these cells, and modulated fractionations with a lower dose of radiation corre- 
sponding to higher radiosensitivity could lead to a reduction of harmful effects. 

With MBED, this study was not intended to implement a finely tuned model based 
on accurate data obtained from preclinical analysis. The aim was to demonstrate the 
potential of the model and its malleability in terms of including further information 
that selective preclinical studies may provide [19]. 

In addition, previous analyses have depended on the validity of the LQ model, which 
has limitations. In particular, the LQ model used during this study does not include the 
time factor. In the generalized LQ model [5,10] the temporal factor is affected by differ- 
ences in EGFr expression due to its influence on potential doubling time, [29-32]. 
This temporal factor can be particularly important when the MBED model is used to 
compare treatment schedules that differ in terms of overall treatment times, tumour 
control or acute effects (where time dependent repopulation may be important). The 
difference of doubling time between the High EGFr group and the Low EGFr group 
identified during the current study will be investigated further in new studies. This dif- 
ference in terms of Tj) can be transformed into an equivalent dose that would be 
required to offset the modified proliferation occurring in one day. The value of this 
equivalent dose can be taken into account during the previous analysis. 
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Overall, in practical applications of the MBED concept, there should be careful con- 
sideration of the relevant physical dose variations, the possible range of biological para- 
meters and pertinent clinical factors. The prudent clinical oncologist should use MBED 
as a guide during clinical decisions rather than as an absolute indicator. The advice of 
acknowledged experts in radiobiological modelling should be sought in more compli- 
cated clinical situations. 

Despite these limitations, the MBED model provides a valid means of accounting for 
modulated intrinsic radiosensitivity effects, which is preferable to neglecting them by 
using a biologically uncorrected physical dose. 

Furthermore, the method is not intrinsically associated with the disease, and can be 
applied to any case by integrating traditional treatment plans and improving the overall 
radiotherapy performances in combined treatments comprising radiosensitizing drugs. 

Conclusion 

During this computational study, the MBED method was introduced. The MBED pro- 
vides a new tool to estimate the effects of heterogeneity on tumour radiosensitivity and 
to assess the dose per fraction required for increased tumour radiosensitivity due to 
EGFr over-expression. Where radiotherapy treatment is combined with radiosensitizing 
drugs, MBED suggests that the fraction sizes modulated according to drug pharmaco- 
kinetics will allow new schedules of dose fractionation to be more effective. 

In conclusion, the MBED method could improve overall radiotherapy performances 
and be utilised to perform more appropriate radiobiological analysis, particularly when 
combined treatment comprising radiation and biological drugs is employed. 

Competing interests 

The authors declare they have r^o competing interests. 
Authors' contributions 

PP developed the model and designed the study. BAJ, LS, BV, DA, MC, FB, Gl checked the appropriateness of the study 
from oncology, radiotherapy and mathematical points of view. PP, RC, MC and LS compiled the manuscript and 
produced the graphical illustrations. AN, GS, MB and RO supervised the manuscript from radiobiological and clinical 
point of view. All co-authors approved the manuscript. 

Acl<nowledgements 

This work resulted from collaborative research between the CROB of Rionero in Vulture and lEO of Milan. CROB 
funded this collaboration. 

We thank Cossu Rocca Maria - Senior Medical oncologist lEO Milan - for expert assistance concerning 
pharmacokinetics. 

Author details 

^I.R.C.C.S. C.R.O.B Regional Cancer Hospital, Rionero in Vulture, Italy. ^U.O. of Radiotherapy, I.E.O. European Institute of 
Oncology, Milan, Italy. ^University of Milan, Milan, Italy. ^Laboratory of Medical Physics and Expert Systems, Regina 
Elena National Cancer Institute, Rome, Italy. ^U.O. of Radiobiology, C.N.A.O, Pavia, Italy. ^Service of Medical Physics, I.E.O. 
European Institute of Oncology, Milan, Italy. ^Service of Medical Physics, Scientific Institute of Tumours of Romagna I.R. 
S.T, Meldola, Italy. ^Department of Radiation Oncology, IRCCS CROB, 1 Padre Pio Street, 85028, Rionero in Vulture, PZ, 
Italy. 

Received: 4 March 2012 Accepted: 26 May 2012 
Published: 19 June 2012 

References 

1. Lammering G: Molecular predictor and promising target: will EGFr now become a star in radiotherapy? 

Radioth Oncol 2005, 74:89-91. 

2. Harari P, Huang SM: Modulation of molecular targets to enhance radiation. Clin Cancer Res 2000, 6:323-335. 

3. Baumann M, Krause M, Dikomey E, Dittmann K, Dorr W, Kasten Pisula U, Rodemann HP: EGFr targeted anticancer 
drugs in radiotherapy: preclinical evaluation of mechanisms, Radioth Oncol 2007, 83:238-248. 

4. Morgan MA, Parsels LA, Kollar LE, Normolle DP, Maybaum J, Lawrence TS: The combination of epidermal growth 
factor receptor inhibitors with gemcitabine and radiation in pancreatic cancer. Clin Cancer Res 2008, 
14:5142-5149. 



Pedicini et al. Theoretical Biology and Medical Modelling 201 2, 9:23 
http://www.tbionned.conn/content/9/1/23 



Page 12 of 12 



5. Fowler JF: The linear quadratic formula and progress in fractionated radiotherapy. Brit J Rad 1989, 62:679-694. 

6. Niemierko A, Goitein M: Implementation of model for estimating tumour control probability for an 
inhomogeneously irradiated tumour. Radioth Oncol 1993, 29:140-147. 

7. Zaider M, Hanin L: Biologically equivalent dose and long term survival time in radiation treatments. Phys Med 

Biol 2007, 52:6355-6362. 

8. Dubray BM, Thames HD: The clinical significance of ratios of radiobiological parameters. Int J Radial Oncol Biol 
Phys 1996, 35:1099-1111. 

9. Barendsen GW: Dose fractionation, dose rate and isoeffect relationships for normal tissue responses. Int J 

Radiat Oncol Biol Phys 1 982, 8:779-790. 

10. Hendry JH, Bentzen SM, Dale RG, Fowler JF, Wheldon TE, Jones B, Munro AJ, Slevin NJ, Robertson AG: A modelled 
comparison of the effects of using different ways to compensate for missed treatment days in radiotherapy. 

C//n Onco/ 1996, 8:297-307. 

11. Zagars GK, Schultheiss TE, Peters LJ: Inter tumour heterogeneity and radiation dose control curves. Radioth 
Onco/ 1987, 8:353-362. 

12. Eriksen JG, Steiniche T, Overgaard J, Danish Head and Neck Cancer Study group: The role of epidermal growth 
factor receptor and E-cadherin for the outcome of reduction in the overall treatment time of radiotherapy of 
supraglottic larynx squamous cell carcinoma. Acta Oncol 2005, 44:50-58. 

13. Bengt KL, Brahme A: The radiation response of heterogeneous tumours. Phys Med 2007, 23:91-99. 

14. Webb S, Nahum AE: A model for calculating tumour control probability in radiotherapy including the effects 
of inhomogeneous distribution of dose and clonogenic cell density. Phys Med Biol 1993, 38:653-666. 

15. Ang KK, Berkey BA, Tu X, Zhang HZ, Katz R, Hammond EH, Fu KK, Milas L: Impact of Epidermal Growth Factor 
Receptor expression on survival and pattern of relapse in patients with advanced head and neck carcinoma. 
Cancer Res 2002, 62:7350-7356. 

16. Dittmann K, Mayer C, Rodemann HP: Inhibition of radiation induced EGFr nuclear import by C225 (Cetuximab) 
suppresses DNA PK activity. Radioth Oncol 2005, 76:157-161. 

17. Akashi Y, Okamoto I, Iwasa T, Yoshida T, Hatashita Y, Yamada Y, Satoh T, Fukuoka M, Ono K, Nakagawa K: 
Enhancement of the antitumour activity of ionising radiation by nimotuzumab, a humanised monoclonal 
antibody to the epidermal growth factor receptor, in non-small cell lung cancer cell lines of differing 
epidermal growth factor receptor status. Br J Cancer 2008, 98:749-755. 

18. Bonner JA, Maihie NJ, Folven BR, Christianson TJ, Spain K: The interaction of epidermal growth factor and 
radiation in human head and neck squamous cell carcinoma cell lines with vastly different radiosensitivities. 
Int J Radiat Oncol Biol Phys 1 994, 29:243-247. 

19. Barendsen GW: Differences in radiosensitivity among cells in culture and in experimental tumours: 
significance for the effectiveness of human cancer therapy. Radioth Oncol 1987, 8:285-300. 

20. Fracasso PM, Burris H III, Arquette MA, Govindan R, Gao F, Wright LP, Goodner SA, Greco FA, Ones SF, Willcut N, 
Chodkiewcz C, Pathak A, Springett GM, Simon GR, Sullivan DM, Marcelpoil R, Mayfield SD, Mauro D, Garrett CR: A 
Phase 1 Escalating Single Dose and Weekly Fixed Dose Study of Cetuximab: Pharmacokinetic and 
Pharmacodynamic Rationale for Dosing. Clin Cancer Res 2007, 13:986-993. 

21. Joiner MC, van der Kogel AJ: The linear-quadratic approach to fractionation and calculation of isoeffect 
relationships. In Basic clinical radiobiology. Edited by Steel GG. London: Arnold; 1997:107-122. 

22. Barendsen GW: RBE as a function of dose for effects on tissues and tumours assessed by the linear quadratic 
model. Int J Radiat Oncol Biol Phys 2000, 46:684-685. 

23. Carlson DJ, Stewart RD, Li XA, Jennings K, Wang JZ, Guerrero M: Comparison of in vitro and in vivo a/(3 ratios 
for prostate cancer. Phys Med Biol 2004, 49:4477-4491. 

24. Begg AC, Vens C: Genetic manipulation of radiosensitivity. Int J Radiat Oncol Biol Phys 2001 , 49:367-371 . 

25. Brown JM: Therapeutic targets in radiotherapy. Int J Radiat Oncol Biol Phys 2001 , 49:31 9-326. 

26. Strigari L, D'Andrea M, Abate A, Benassi M: A heterogeneous dose distribution in simultaneous integrated boost: 
the role of the clonogenic cell density on the tumour control probability. Phys Med Biol 2008, 53:5257-5273. 

27. Roberts SA, Hendry JH: A realistic closed-form radiobiological model of clinical tumour-control data 
incorporating intertumour heterogeneity. Int J Radiat Oncol Biol Phys 1998, 41:689-699. 

28. Thariat J, Milas L, Ang KK: Integrating radiotherapy with epidermal growth factor receptor antagonists and other 
molecular therapeutics for the treatment of the head and neck cancer. Int J Radiat Oncol Biol Phys 2007, 69:974-984. 

29. Suwinski R, Jaworska M, Nikiel B, Grzegorz W, Bankowska-Wwozniak M, Wojciech M, Krysztof S, Dariusz L: 
Predicting the effect of accelerated fractionation in postoperative radiotherapy for head and neck cancer 
based on molecular marker profiles: data from a randomized clinical trial. Int J Radiat Oncol Biol Phys 2010, 
77:438-446. 

30. Chung CH, Zhang Q, Hammond EM, Trotti AM III, Wang H, Spencer S, Zhang HZ, Cooper J, Jordan R, Rotman MH, 
Ang KK: Integrating epidermal growth factor receptor assay with clinical parameters improves risk 
classification for relapse and survival in head-and-neck squamous cell carcinoma. Int J Radiat Oncol Biol Phys 

2010,81:331-338. 

31 . Bentzen SM, Atasoy BM, Daley FM, Dische S, Richman PI, Saunders Ml, Trott KR, Wilson GD: Epidermal growth factor 
receptor expression in pretreatment biopsies from head and neck squamous cell carcinoma as a predictive factor 
for a benefit from accelerated radiation therapy in a randomized controlled trial. J Clin Oncol 2005, 23:5560-5567. 

32. Smid EJ, Stoter TR, Bloemena E, Lafleur MV, Leemans CR, van der Waal I, Slotman BJ, Langendijk JA: The 
importance of immunohistochemical expression of EGFr in squamous cell carcinoma of the oral cavity 
treated with surgery and postoperative radiotherapy. Int J Radiat Oncol Biol Phys 2006, 65:1323-1329. 



doi:l 0.1 1 86/1 742-4682-9-23 

Cite this article as: Pedicini et al.: Modelling the correlation between EGFr expression and tumour cell 
radiosensitivity, and combined treatments of radiation and monoclonal antibody EGFr inhibitors. Theoretical 
Biology and Medical Modelling 201 2 9:23. 



